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ABSTRACT 

We  discuss  the  flow  of  BKZ  fluids  in  an  orthogonal  rheometer.  Some  analytical  results 
are  proved,  and  numerical  solutions  are  obtained  for  the  Currie  model.  These  solutions 
show  a  boundary  layer  behavior  at  high  Reynolds  numbers  and  the  possibility  of  discon¬ 
tinuous  solutions  or  nonexistence  at  high  Weissenberg  numbers.  . 
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SIGNIFICANCE  AND  EXPLANATION 


The  orthogonal  rheometer  is  a  device  used  in  rheological  measurements  which  con¬ 
sists  of  two  parallel  plates  rotating  with  the  same  speed  about  parallel  distinct  axes.  In 
this  paper,  the  equations  of  motion  for  flows  in  such  a  rheometer  are  investigated,  with 
particular  emphasis  on  the  possibility  of  loss  of  uniqueness.  A  particular  model  is  stud¬ 
ied  numerically,  and  nonunique  solutions  with  discontinuous  derivatives  are  found.  The 
relevance  of  such  a  behavior  in  relation  to  experiments  is  discussed. 
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FLOW  OF  VISCOELASTIC  FLUIDS  BETWEEN  PLATES 


ROTATING  ABOUT  DISTINCT  AXES 


K.R.  RajagopaP,  M.  Renardy,  Y.  Renardy  and  A.S.  Wineman^ 


1.  Introduction 

The  flow  in  an  orthogonal  rheometer  has  been  the  object  of  considerable  interest  and 
has  been  studied  by  several  authors  (cf.  Maxwell  and  Chartoff  [l],  Blyler  and  Kurtz  \2\, 
Bird  and  Harris  l3  ,  Huilgol  4  ,  Kearsley  5  ,  Gordon  and  Schowalter  [6]).  The  apparatus 
consists  of  two  parallel  plates  which  rotate  with  the  same  constant  angular  velocity  about 
distinct  axes  (see  Fig.  1).  In  all  these  investigations  the  analysis  was  based  on  the  assump¬ 
tion  that  inertial  effects  are  negligible.  Abbot  and  Walters  j7;  obtained  an  exact  solution 
for  the  flow  of  a  Newtonian  fluid  under  the  assumption  that  the  plates  are  infinite.  They 
also  studied  the  flow  of  a  viscoelastic  fluid  in  such  an  apparatus,  but  in  this  case  they 
assumed  that  the  distance  between  the  two  axes  of  rotation  is  small  and  employed  a  power 
series  expansion  in  the  distance  between  the  axes.  This  kind  of  expansion  is  also  studied 
by  Goldstein  and  Schowalter  8;. 

Rajagopal  9  showed  that  the  equations  of  motion  for  an  arbitrary  simple  fluid  are 
compatible  with  the  following  ansatz  for  the  velocity  field 

u  ^  -  n'y  ■  g(z)], 

v^ax-f(z)\  (1) 

u’  =  0. 

Here  u.r  and  u'  are  the  x.y  and  z  components  of  the  velocity,  and  the  plates  of  the 
rheometer  are  parallel  to  the  j-,y-plane.  This  velocity  field  describes  a  flow  in  which  each 
plane  rotates  rigidly  with  the  same  angular  velocity  fi.  The  locus  of  the  centers  of  rotation 
is  in  general  a  curve  in  space  defined  by  x  f(z)-  y  3(2)-  The  flow  given  by  (1)  is 
a  motion  with  constant  stretch  history  (cf.  Coleman  10.  Noll  !lT).  This  implies  that 
the  equations  of  motion  for  a  simple  fluid  reduce  to  the  same  order  as  the  Navier-Stokes 
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equations  and  that  the  no  slip  boundary  conditions  are  sufficient  for  determinacy.  For 
related  results,  see  also  Goddard  [12|. 

For  certain  special  fluid  models,  the  resulting  boundary  value  problem  can  be  solved 
exactly.  Rajagopal  and  Gupta  [13]  obtained  an  exact  solution  for  the  incompressible  fluid 
of  second  grade,  while  Rajagopal  and  Wineman  [14:  established  exact  solutions  for  a  special 
class  of  BKZ  fluids. 

In  this  paper,  we  study  the  flow  of  a  general  BKZ  fluid  in  the  orthogonal  rheometer. 
Detailed  numerical  results  are  provided  for  a  specific  model.  We  find  that  it  is  possible 
that  the  locus  of  the  centers  of  rotation  is  piecewise  continuously  differentiable. 

2.  Equations  of  motion 

The  Cauchy  stress  T  in  an  incompressible  homogeneous  simple  fluid  undergoing  a 
motion  with  constant  stretch  history  can  be  expressed  as  (cf.  Wang  [15]) 

T  = -pl  +  f(A,,A2,A3),  (2) 

where  A].  A2  and  A3  denote  the  first  three  Rivlin-Ericksen  tensors  (cf.  Rivlin  and  Erick- 
sen  [16  )  defined  recursively  by 

A,  =  (Vu)  + (3)1 


A„  =  ^^  4  A„-,(Vr)  +  (Vu)^A„_,.  (3)2 

In  these  equations,  v  denotes  the  velocity  vector,  ^  is  the  material  time  derivative  and 
Vv  is  defined  according  to  the  convention  dv^jdxj. 

For  the  motion  under  consideration,  it  can  be  shown  that  (cf.  Rajagopal  [9]) 

A3  -  n'A,.  (4) 

and  hence  (2)  reduces  to' 

T-  pl^f(A,.A2).  (5) 

Since  Aj  and  A 2  depend  only  on  z,  the  equation  of  motion  reduces  to 

d/ 1 3  -  d/ 2,3  T  d/ 33  ■ 

—  i  -  }  ‘  .  k 

dz  dz  dz 

Vp  .  V(p<I>)  pD*  T  f{z)  i  y  gizYl  (6) 

The  same  form  is  obtained  for  v'iscornetric  flows  where  A3  0.  However,  the  function 
f  ap[)earing  in  (5)  is  in  general  different  from  that  in  viscometric  flow. 
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where  i,  j  and  k  denote  unit  vectors  in  the  i,  y  and  z  directions,  respectively.  In  obtaining 
(6)  we  have  also  assumed  that  the  body  force  is  conservative  and  derivable  from  a  potential 
By  operating  on  (6)  with  the  curl  operator,  we  find 


(7)i 

'J"  = 

(7)2 

These  equations  can  be  integrated  to  yield 

(8)1 

~  9, 

dz 

(8)2 

where  q  and  s  are  constants.  If  /  and  g  satisfy  equations  (8),  then  (6) 
only  if  the  pressure  is  given  by 

is  satisfied  if  and 

P=  -  +  ^(x^  +  y^)  ^  sx  +  qy  ^  733(2)  -t-  C, 

(9) 

where  C  is  an  arbitrary  constant.  Hence  the  system  (8)  can  be  regarded  as  equivalent  to 

(6). 

Since  the  fluid  adheres  to  the  rotating  plates,  the  appropriate  boundary  conditions 
for  the  problem  are  (see  Fig.  1) 

w  =  —  -  fly,  V  =  fix,  w  =  0  at  z  =  h,  (10)  i 

2 

u  ='  -  —  -  riy,  V  =  fix,  w  -  0  at  z  =  0.  (10)2 

2 

When  inserted  into  (1),  these  conditions  lead  to 

/(o)  =  /(fi)  =  0,  g(0)  -  9(/i)  =  (11) 

For  the  following,  we  consider  fluids  in  which  the  Cauchy  stress  tensor  T  is  related  to 
the  history  of  the  deformation  in  the  following  manner  (cf.  Bernstein,  Kearsley  and  Zapas 
17) 

T  ^  -pi  -  2^  {r,C,-'(r)  ~  f'2Ct(r)}  dr.  (12) 

Here  C,  is  the  right  relative  stretch  tensor  defined  by 

C,(7)  ^ 
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(13) 


and  Ut  is  the  derivative  of  the  stored  energy  function  U  =  -  t 

the  argument  The  latter  are  defined  by 

)  with  respect  to 

/,  =  trCr'(r),  /2  =  tr  Q(r). 

(14) 

It  follows  from  a  lengthy  but  straightforward  computation  that  the  equation  of  motion 
reduces  to 

J 

^{/'(2)fi(K)  +  g[z)A{K)}  pCl^fiz)  +  s. 

(I5)i 

-  f'(z)A(K)  4  9'(2)B(/c)}  pn^g(z)  +  q. 

(15)2 

The  functions  A  and  B  are  defined  by 

r  OO 

A(k)  ~  2  1  f/[3  +  2(l  -  cosflo)«:,  oi  sin  flo  da, 

Jo 

(I6)i 

B{k)  2  f  ^'[3  +  2(1  -  cos  na)«:,  al(l  -  cos  flo)  do. 
Jo 

(16)2 

where 

U J,a]  =  ih{lJ,a)  U2(IJ,a). 

(17) 

The  quantity  k  is  defined  as 

K{z]  r.  f{zy^  -  [g'{zy^. 

(18) 

In  this  work,  we  shall  find  numerical  solutions  of  the  system  (15)  with  boundary  conditions 
(11)  for  a  special  class  of  BKZ  fluids.  Before  we  do  this,  however,  we  first  remark  on  a  few 
general  characteristics  of  the  problem  under  consideration. 


3.  Analysis 


The  solution  to  the  system  (15)  depends  on  the  material  properties  of  the  model. 
However,  when  a  q  s  0.  we  always  have  the  trivial  solution  /  -^  g  -  0.  If  a,  q 
and  .s  are  small,  this  trivial  solution  can  be  uniquely  continued  to  a  solution  f(z\a,q,s) 
and  g{z:a,q.s).  .4s  the  values  of  a.  q  and  s  are  further  increased,  the  process  of  unique 
continuation  can  only  fail  if  one  of  the  following  situations  occurs: 

(i)  th(‘  equations  become  singular. 

(ii)  the  solution  and  or  its  derivative  become  infinite. 

(iii)  the  linearization  of  (15)  at  f{z\a.q.s).  g(z:a.q,s)  has  a  nontrivial  solution. 

VVhethc'r  any  of  these  situations  obtains  depends  on  the  choice  of  the  constitutive  equation 
wliich  determines  tfie  functions  ,1  and  li.  In  the  following,  we  discuss  some  implications 
of  each  of  these  conditions. 
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Equations  (15)  imply  that 

r\B  +  2B'f'^  +  2A'f'g'\  +  g"\A  +  2B'f'g'  -  2A'g'^\  =  pH V  ^  (19)  i 

r\-A  -  2A7'^  -  2B7'9'i  +  -  2A' f'g' ^  2B'g'^\  =  +  9.  (19)2 

This  system  becomes  singular  if  and  only  if  the  determinant  of  the  coefficients  of  the  second 
derivatives  becomes  zero.  This  determinant  is  given  by 

A^-  B^^  2k{AA'  -  BB')  =  ^{k(A^  -  B^)}.  (20) 

aK 

If  the  function  k  —  k(j4^  B^)  is  strictly  monotone,  then  D  can  never  be  zero.  This  is 

the  case  for  the  models  studied  by  Rajagopal  and  Gupta  J3i  and  Rajagopal  and  Wineman 
14  where  A  and  B  are  constant. 

Case  (it) 

We  now  assume  that  the  system  (15)  is  non-singular,  and  discuss  some  of  the  implica¬ 
tions  of  (ii).  The  following  lemma  shows  that  if  D  and  B  are  positive,  and  the  derivative 
of  /  or  9  becomes  infinite,  then  this  must  happen  at  the  boundary. 

Lemma:  If  D  and  B  are  positive,  then  k  cannot  have  an  interior  maximum. 

Proof:  If  we  differentiate  equations  (15)  and  evaluate  at  k'  —  0,  we  find  that 

/"'B  *  9'"^  +  /'BV"  ^  g'A'x'*  =  (21), 

-r'A  -  g"'B  -  f'A'K"  +  g'BV"  =  pQ^g’.  (21)2 

Multiplying  (21),  by  f'B  -  g'A  and  (21)2  by  g'B  -  f'A,  and  adding  the  two,  we  obtain 

(A^  ^  B^){f"'f'  +  g"'g')  -  {AA'  ^  BB’)kk"  -  pQ^B.  (22) 

It  follows  from  the  definition  (18)  of  k  and  from  (20)  that  (22)  can  be  rewritten  as 

^-Dk"  =  {A^  -  -  g"'')  -  pn\B.  (23) 

Since  B  and  D  are  positive,  it  follows  that  k"  >  0.  Hence  any  interior  point  where  fc'  =  0 
is  a  minimum. 

If  further  assumptions  are  added,  it  is  possible  to  say  more. 

Lemma:  If  A  and  B  are  positive,  bounded  functions  of  n  and  they  have  strictly 
positive  lower  bounds,  then  case  (ii)  does  not  happen. 

Proof:  W’e  multiply  (15),  by  /,  (15)2  by  g,  add  them  and  integrate.  This  leads  to 

j  pQ'if^  -  5/7  "  ^(/'‘  *  y'^) '  -7  *  yy  ' 


>Vy‘J 


•I!"  .".  A  /• 


l[-f'{h)A(K{h))  g'ih)B{K{h))  -  rmiKiO))  -  g'{0)B{K{0)),  (24) 

Next,  we  multiply  (15) i  by  g  and  (15)2  by  /,  subtract  them  and  integrate.  This  leads  to 


^  g'^)  *  sg  -  qf  dz  ^ 


^\f'(h)B{K(h))  *  g'(/i),4(/c(/i))  ^  r(0)B(K(0))  *  ff'(0).4(/c(0)):.  (25) 

It  follows  easily  that  the  //‘-norm  of  /  and  g  can  be  estimated  in  terms  of  the  boundary 
terms  occurring  on  the  right  side  of  either  (24)  or  (25).  If  the  //‘-norm  stays  bounded 
while  following  a  branch  of  solutions,  then  the  minima  of  f'B  •  g'A  and  -  f'A  g'B 
also  stay  bounded,  and  so  do  the  maxima  of  /  and  g  .  By  integrating  (15)  with  respect 
to  z.  it  follows  that  f'B  -  g'A  and  f'A  -  g' B  remain  uniformly  bounded,  and  the 
assumptions  on  A  and  H  imply  that  /'  and  g'  also  remain  uniformly  bounded.  In  order 
for  case  (ii)  to  happen,  it  is  therefore  necessary  that  the  terms  on  the  right  of  (24)  and 
(25)  both  tend  to  infinity.  Again  by  integrating  (15)  with  respect  to  2,  we  find  that  this  is 
possible  only  if  the  minima  of  f  B  ■  g'A  and  f'A  -  g'B  both  tend  to  ••  oc.  In  this  case, 
however,  we  conclude  that 


tends  to  •  oc  uniformly  on  O.h  .  which  is  inconsistent  with  the  boundary  conditions. 

Whether  case  (ii)  ever  obtains  depends  of  <ourse  on  the  model  under  consideration. 
In  the  models  considered  by  Rajagopal  and  (lupta  1.3  and  Rajagopal  and  Wineman  14  , 
the  above  lemma  shows  that  case  (ii)  cannot  obtain. 

Case  (iii) 

Depending  on  the  model  under  consideration,  this  case  could  lead  to  limit  points  or 
bifurcations.  This  cannot  happen  as  long  as  the  matrix  of  the  coefficients  of  f"  and  g"  in 
equations  (19)  remains  positive  definite.  This  is  the  case  as  long  as 

B  '  0  and  //‘  -  2BB'k  A''\'^  >  0.  (27) 


If  A  and  B  are  positive  constants,  such  as  in  13  and  14  ,  (27)  holds." 

In  summary,  we  found  that  neither  of  the  cases  (i)-(iii)  ever  occurs  in  a  simple  model 
where  A  and  B  are  positive  constants.  Below,  we  report  numerical  results  for  a  model 
where  A  and  B  are  rather  complicated  expressions.  The  results  clearly  indicate  that  case 
(i)  occurs  in  this  model. 

■  Of  course  this  analysis  only  shows  the  absence  of  bifurcations  within  the  class  of 
solutions  described  by  the  kinematics  of  equation  (1).  rh«Te  itia\  well  be  other  possible 
flows  which  cannot  be  desribed  by  (I). 


(i 
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-  • 


4.  Numerical  Solutions 


Let  the  interval  [0,h\  be  partitioned  into  N  equal  subintervals  bounded  by  points 
2.  =  (i  -  l)h/N,  i  ^  1,2,. ..,N  +  1.  Define  the  notation  /,  =  /{z,),  9i  ^  g{zi)-  The  finite 
difference  scheme  used  to  approximate  the  system  (15)  is  given  by 


A"(A-/  B(k-)  +  A-g  >!(«-))  + 

(28), 

A'"  (- A~/  v4(/c_)  ^  A~g  B(#c-))  =  +  q. 

(28)2 

where,  with  Ah  —  h'S, 

(29) 

and 

-  (A-f)^-(A-9)\ 

(30) 

Since  by  (11),  /,  =  =  0  and  3,  -  -a/2,  =  a/2,  there  are  2N  -  2  unknowns, 

/j,  fs  and  g2 .  gs-  This  discretization  leads  to  the  following  nonlinear  equations 


(ft~i  -  (g.^i  -  g,)A(K.^)  -  (/,  -  /,-i)B(k.  ' 

-(g,  -  g,.  ,)>!(«-)  -  (A/i)2(pDV.  -  -  0,  (31), 

-  (/.-i  -  /.)^(«+)  (9.-1  -  gt)B{K+)  ••  (/,  -  /,_,)>4(«_) 

-[g,  -  p,_,)B(k_)  -  (Ah)^pn^g,  -  9)  =  0  (31)^ 

for  each  mesh  interval  point  2,,  i  =  2,  ...,  N. 

Denote  equations  (31),  and  (31)2  respectively  by 

E\(fi-i,g,-  1,  /j, 9i, /i- ii9»+i)  =  e[  ^  -  0,  (32) , 

B2(/i  -  1, 9i  -  i>  /i< /«-  !•  9i  + 1)  =  Bj  *  =0,  (32) 2 

i  -  2,  N. 

This  system  is  solved  using  Newton's  scheme  on  the  vector  equation  V'  =  0  where 


The  corresponding  Jacobian  is  {2N  -  2)  square  and  is  banded  with  three  sub-diagonals 
above  and  below  the  main  diagonal: 


dB\-' 

dEl'' 

dE[-' 

0 

df. 

dg^ 

dfz 

dgz 

c3 

ag/-' 

dEl"' 

dE^-' 

0 

dh 

dg2 

dfz 

dgz 

at-;/' 

dE'/' 

dE'/' 

dE\'' 

dE'^' 

3f, 

dg. 

3h 

dg-. 

dg* 

ai,'!'’ 

dE^-' 

dE'S' 

a  El'-’ 

dEl'' 

dE\y 

dg^ 

3fz 

dgz 

dU 

dg* 

0  0 
0  0 


0 

I  0 


0 

0 


dE'J^' 

dEl"' 

dE\'^' 

37.V  I 

dgs  -  1 

dfN 

dgs 

dE\^  ' 

dEi''’’ 

dE[^' 

dE!,'^' 

a/N-T 

1 

i  k. 

dfN 

dgs 

(34) 


Several  tests  were  carried  out  for  the  convergence  of  Newton’s  iteration  scheme.  The 
first  uses  the  exact  solution  obtained  by  Rajagopal  and  Wineman  [14]  for  the  czise  in  which 
.‘\(k)  and  /f(K)  are  constants  and  s  =  q  —  0.  Let  and  denote  the  values  of  /(z) 
and  g(z)  given  by  the  exact  solution  in  il4i  and  let  denote  the  values  computed 

from  .Newton's  iteration.  Calculations  were  carried  out  for  p  -  3.  fl  ~  1,  a  =  0.1,  h  = 
2,  3,  -  2  and  a  .Newton's  scheme  tolerance  of  10  The  absolute  errors 

’  /,'  '•  9^'  s!  '  t'valuated  as  the  number  of  points  was  doubled  from  N  =  10  to 

.V  -  20  and  doubled  again  to  .V  40.  The  errors  went  down  by  1/4  each  time  indicating 
quadratic  convergence  as  would  theoretically  be  expected. 

Next.  A(k)  and  B(k)  were  prescribed  to  be 


.4  ( K  )  o  I  K  -  o  2  • 


(35), 


//(a)  •  62.  (35)2 

with  n  hi .  h-j  prescribt'd  cotistants.  .\n  exact  solution  for  /(c)  and  g(c)  is  unknown  for 
thi^-  case.  In  order  to  test  the  converg<'nce  f)f  tlie  Newton  iteration  scheme,  the  problem 
was  modified  in  the  following  manner.  /(;)  and  9(c)  were  chosen  in  the  form 


/(^)  /,((-'  h  -ir  \h  2)-').  (36), 

y(-)  </i(-  h  '-'I'  v.i-  h  2).  (36)2 

where /,.(/ 1 art' prescribed  (oust  ant  s.  Note  that  in  at  (  ordam  e  wit  h  (II). /(O)  f(h] 

0  ant]  f/(())  q{h).  The  /  ami  g  as  gi\en  b\  (.36i  do  not  saii^ft  t'tjuations  (15).  W'e  ran. 

tif>wi'v<'r.  make  tht'iri  satisfy  the  equatit)ns  b\  adding  a  ftirt  mg  itTrn  to  the  right  hand  sitle 
of  (15).  \\  ht  II  tlist  ret izing.  we  obtain  a  t orre-ptinding  ftirt  mg  term  tti  be  atitietl  tt)  the 

right  side  of  (32).  We  dentite  these  terms  b\  /, ,  .  . 


The  original  system  (32)  is  thus  modified  to 


=  E\ 


(«) 


(37)  1 


£2"  “  '  2 . A’.  (37)2 

for  which  (36)  1  and  (36)2  provide  the  exact  solution.  Computations  were  carried  out 
for  the  two  sets  of  parameters;  (p,n,/i,ai,a2'^i-^2-/iiffii92)  —  (3. 1,2, 1,1, 1,1, 1,1,1). 
(3,  1,  2,  1.5.  2.9.  3.1.  10.0,  0.1.  2.4,  5.2)  with  10.  20,  40  and  80  mesh  points.  Comparison 
of  the  computed  values  and  its  exact  solutions  showed  quadratic  convergence. 

In  the  final  case,  the  process  was  repeated  for 


A(k) 

2 

0  1 K  •  flo  *  ^3^  » 

(38) 

B(k] 

6ik  •  62  * 

(38) 

( 'omput  at  ions  were  performed  for  the  second  choice  of  parameters  above,  with  03  =  1,  63  - 
0  and  with  a..  1,  f),-,  1.  The  results  again  showed  quadratic  convergence. 


5.  The_Currie_  .Model 

Computations  were  carried  out  using  the  form  of  the  stored  energy  function  V  derived 
by  Currie  18  as  an  approximation  for  the  Doi-Edwards  model  (cf.  19  ).  It  is  given  by 

^  rn(s)0{J),  (39) 

where 

V{J)  5ln(J  -  1)  -  9.73. 

J  ^  h  ^  2(/2  '  3.25)'^^ 

and 

m(.s)  f7(.s). 

t7(.s)  is  the  stress  relaxation  function  of  linear  viscoelasticity 

m(^)  ~  C'e' 

The  form  of  the  function  f’(/.a)  defined  in  (17)  is 

C(I.a]  re  ■''*/•(/)  (44) 

where 


(40) 

(41) 

(42) 

.  For  the  present  calculations, 

(43) 


r{i] 


,> 

/  1  •  2\  /  .3.25 


(45) 


9 


Functions  A  and  B  defined  through  (16)  become 


A  =  2C  I  +  2(1  -  cosna)/c)  sin  ria  da, 

Jo 

B  =  2C  f  2(1  -  cosno)K)(l  -  cos  fla)  da. 

J  (.1 


(46), 


(46): 


For  convenience  in  numerical  integrations,  these  integrals  can  be  re-expressed  as 


2C 


B  = 


n(l  -  exp(-27rA/n)) 
2C 


n(l  -  exp(-  27rA  Tl)) 


L 


L 


2ir 


2n 


e  “L' (3  -  2(1  -  cos  a*)K:)  sin  a' da^ 


Xck*  /n f : 


t/ (3  4  2(1  -cosq')k)(1  -  cos  o')  do'. 


(47). 

(47)2 


where  the  non-dimensional  time  a'  =  off  has  been  introduced. 

The  parameters  of  the  problem  are  now:  a,h  with  dimension  L;  12.  A  with  dimension 
T~^;  p  with  dimension  \4 1  L^\  and  C  with  dimension  M:[LT^).  Note  C/A^  is  the  zero 
shear  rate  viscosity.  From  these  the  following  non-dimensional  parameters  are  formed; 
alh,  a  Weissenberg  number  H  e  -  12  A  and  a  Reynolds  number  Re  ph^ilX^ jC .  These 
appear  in  the  non-dimensionalized  version  of  equations  (15) 


'  (b%  -  )  -  Re  H  e  /  •  .s. 

z  \  dz  dz  ' 


d 

dz 


ReWeq^q, 

\  dz  dz  ' 


(48). 

(48)2 


where  /  -  //a,  g  g/a,  z  =  zjh.A  -  /lA  C  and  B  BXiC. 

The  Jacobian  eissociated  with  Newton's  method  requires  the  evaluation  of  the  deriva- 
tives  ^4  and  .  Expressions  for  these  are  given  by 

d  A  A  C'*  r  2 

-  J  e  “f’'(.3  *  2(1  coso')s')(l  -  cos  o')  sin  o’ do^ 


dK  12(1  exp(  27rA  12)) 

dB  1C 

dx  12(1  exp(  27.  A  12)) 


(49), 


/ 


-  Hr'' 


L’’(3  *  2(1  coso’)k)(1  coso’)^  da'.  (49)2 


Each  of  the  integrals  in  (47)  and  (49)  is  evaluated  by  subdividing  the  interval  [0.27ri 
into  subintervals  of  equal  length  and  then  using  the  trapezoidal  rule. 

The  following  are  the  components  of  traction  at  the  walls  z  O.h  which  were  also 
computed: 

C  B(x.n)f'  ■  .■\(x.U)g'.  (50), 


.4(7c.12)/'  •  B{x.n)g'. 


(50)2 


C-;': 

'S'l* 


*  •. 


•  *  . 

w  *  ■ 

■  •.  *, 

r: 
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A  ^ 


By  comparing  this  with  equations  (15),  we  find  that  if  sh  +  pfl*  /(zjdz  is  not  zero,  then 
the  component  of  traction  tx  has  different  values  at  the  top  and  bottom  plates.  Similarly, 
if  qh  +  f^g{z)dz  is  not  zero,  then  the  component  of  traction  ty  would  be  different  at 
the  top  and  bottom  plates.  This  was  also  checked  numerically  for  the  Figures  discussed 
below. 


6.  Results  and  Discussion 

The  three  non-dimensional  parameters  governing  the  problem  are  the  Reynolds  and 
the  Weissenberg  numbers  and  a/h.  Rajagopal  and  Wineman  |14]  restricted  themselves  to 
the  special  case  when  s  =  q  =  0.  In  the  case  of  small  Reynolds  numbers,  they  found  their 
solutions  for  /  and  g  to  be  approximately  linear.  However,  at  high  Reynolds  numbers, 
they  found  a  boundary  layer  structure  in  their  solution.  The  model  considered  here  is 
much  more  complicated  than  that  considered  in  il4l  and  it  is  possible  that  the  solutions 
possess  a  more  complicated  structure. 

We  shall  first  discuss  the  situation  when  the  Reynolds  number  is  small  and  the  Weis¬ 
senberg  number  is  large.  When  s  -  ^  =  0,  we  find  that,  below  some  critical  value  of  a/h, 
the  numerical  scheme  converges  to  a  solution  w'hich  is  anti-symmetric  about  the  mid-plane 
with  /  =s  0  and  g  being  approximately  linear  as  indicated  in  Figures  2  through  5.  (Note 
that  f  -  0  and  g  linear  is  an  exact  solution  for  Re  =  0.)  Figure  2  shows  that  /  ~  fl 
for  small  fl  and  /  ~  a  for  small  a.  Figure  3  shows  that  g(z)  ~  a(z/h  -  0.5)  for  small 
n  and  a.  At  We  ~  10,  Re  -  2.  Newton’s  scheme  converged  quickly  until  a  reached  ap¬ 
proximately  .94,  above  which  it  diverged.  As  the  critical  value  of  a  is  approached,  the 
discriminant  D  =  ^  +  2k{AA'  -r  BB')  approaches  zero  and  the  Newton  scheme 

has  convergence  problems.  The  solutions  shown  have  been  tested  for  convergence  with 
20  and  40  meshpoints.  Figures  4  and  5  display  solutions  with  discontinuous  derivatives 
when  We  -  4,  Re  -  0.4.  The  solutions  possess  the  anti-symmetric  structure  for  a  up  to 
about  0.9750.  above  which  the  discriminant  D  approaches  zero.  For  larger  values  of  a, 
solutions  are  not  unique.  They  may  possess  discontinuous  slopes  with  the  discriminant  D 
becoming  negative  on  part  of  the  solution.  Which  solution  the  Newton  iteration  would 
converge  to  depends  on,  e.g.,  the  number  of  meshpoints.  We  display  the  results  obtained 
with  20  meshpoints  for  a  above  0.9750.  It  is  found  that  the  discontinuous  solutions  are 
not  anti-symmetric.  This  breakdown  in  symmetry  is  purely  due  to  the  round-off  error  in 
the  numerical  scheme  rather  than  any  non-symmetric  input  into  the  problem.  In  view 
of  experiments,  it  is  interesting  to  note  that  the  tractions  on  the  plates  can  be  different 
for  these  discontinuous  solutions.  It  is  evident  from  Figure  4  that  the  integral  of  /  is  not 
zero,  and  hence  tj  must  have  different  values  at  the  two  plates.  Direct  computation  of  the 
tractions  showed  that  both  tj  and  iy  have  different  values  on  the  plates.  Note  that  there 
is  no  difference  in  the  tractions  as  long  as  either  the  solution  is  antisymmetric  or  inertial 
effects  are  ignored. 

When  .s  0  or  9  0.  the  functions  /  and  g  are  not  symmetric,  as  show'n  in  Figures  6 
and  7.  where  .s  -  0.1.  q  0.01.  We  4.  Re  0.4.  Since  our  numerical  scheme  is  designed 
to  start  with  a  -  0.  and  a  is  incremented  in  small  steps,  and  the  solutions  for  /  and  g  are 
obtained  from  iterating  on  the  respective  solutions  for  the  previous  value  of  a,  the  scheme 


r-v-ri 


II 


converges  to  a  smooth  solution  while  a  remains  small.  For  example,  in  Figures  6  and  7,  we 
find  that  the  solutions  are  smooth  until  a  is  slightly  over  0.8.  A  critical  value  of  a  is  then 
reached  at  which  the  discriminant  is  zero,  and  the  scheme  has  convergence  problems.  The 
solutions  plotted  for  a  <  0.8  can  be  reproduced  with  finer  meshes.  Above  the  critical  value 
of  a.  the  Newton  scheme  again  converges  but  the  solutions  are  non-unique  and  can  be 
non-smooth  with  discontinuous  derivatives,  with  the  discriminant  being  negative  on  part 
of  the  solution.  Which  solution  the  .Newton  scheme  chooses  to  converge  to  depends  on  the 
number  of  rrieshpoints  and  may  be  influenced  by  round-off  errors  as  above.  Using  twenty 
meshpoints,  the  scheme  converged  at  a  -  0.9510626  to  the  non-smooth  solutions  displayed 
in  Figures  6  and  7.  The  graph  of  this  solution  has  been  plotted  with  interpolation  from 
the  solution  computed  with  twenty  meshpoints  up  to  z!h  --  0.95.  The  kink  in  the  solution 
occurred  at  the  meshpoint  next  to  the  boundary.  This  is  not  surprising,  since  k  assumes 
its  largest  values  at  the  boundary  (cf.  the  lemma  in  section  3),  and  hence  D  would  first 
become  zero  at  the  boundary. 

The  appearance  of  discontinuous  solutions  may  be  related  to  local  yielding,  which  is 
observed  experimentally  (cf.  Wissbrun  20  ).  We  note  that  although  the  flow  under  study 
is  not  viscometric,  there  is  a  qualitative  similarity  with  simple  shear  flow.  The  Currie 
model  has  a  shear  stress  maximum  in  shear  flow,  and  if  shear  rates  exceed  a  critical  value, 
then  discontinuous  solutions  are  also  possible.  In  the  geometry  studied  in  this  paper, 
solutions  with  discontinuous  derivatives  were  shown  to  exist  for  nonlinear  elastic  materials 
with  non-convex  stored  energy  functions  (cf.  Rajagopal  and  Wineman  21  ). 

If  the  Weissenberg  number  is  small  and  the  Reynolds  number  is  large,  there  is  a 
boundary  layer  structure  to  the  solution  as  is  to  be  expected  from  the  study  of  simpler 
models  13  ,  14  .  see  also  Joseph  22:,  where  the  case  of  large  H  and  small  ajh  is  considered 
for  arbitrary  simple  fluids.  Kxarnples  are  shown  in  figures  8  and  9. 
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Fig.2 

Graph  of  f{zlh)  versus  z/h,  s  =  q  =  0. 

A  =  0.5,  C  =  2.5,  p  =  1,  /i  =  2,  zero  shear  rate  viscosity  =^10. 

(1)  a  =  1.0,  =  0.2,  We  =  1.0. 

(2)  a  =  0.1,  Re  =  2.0,  =  10.0. 

(3)  a  =-  1.0,  --  0.04,1^6  -  0.2. 

(4)  a  -  0.l,Re  =  0.2,  W’e  ^  1.0. 

(5)  a  =  0.l,Re  =  0.04,  Wc  =  0.2. 

Fig.3 

Graph  of  g{z/h)  versus  z  'h,  s  ~  q  -  0. 

A  =  0.5,  C  -  2.5,  p  i,  h  —  2,  zero  shear  rate  viscosity  ==  10. 

_  :  Re  -  0.2,  W’e  -  1.0  and  Re  -  0.04,  IVe  -  0.2. 

. :  Re  -  2.0,We  -  10.0. 

Fig.4 

Graph  of  f{z/h)  versus  z/h,  s  =  q  -  0. 

A-1,  C  ~  1,  p  -  0. 1,  h  =  1,  n-4,  li'e-4.  Re  —  0.4,  zero  shear  rate 

viscosily  =  1.  Twenty  meshpoints  were  used  in  the  finite-difference  scheme. 

Fig.5 

Graph  of  g(z/h)  versus  z/h,  s  =  q  -  0. 

A  --  1,  C  -  1,  p  0.1,  h  —  I,  fl  =  4,  We  =  4,  Re  =  0.4,  zero  shear  rate 

viscosity- 1.  Twenty  meshpoints  were  used  in  the  finite-difference  scheme. 

Fig.6 

Graph  of  f{z:h)  versus  z/h,  s  =  0.1,  q  =  0.01. 

A  =  1,  C  -  I,  p  -  0.1,  h  -  1,  fl  =  4,  We  =  4,  Re  =  0.4,  zero  shear  rale 

viscosity=l. 

Fig. 7 

Graph  of  g^z/h)  versus  z/h,  s  =  0.1,  q  —  0.01. 

A  =  1,  C  -  1,  p  -  0.1,  h  =  1,  n  =  4,  We  -  4,  Re  -  0.4,  zero  shear  rate 

viscosity— 1. 

Fig.8 

Graph  of  f(z/h)  versus  z/h,  .s  --  9  =  0. 

A  -  0.5,  C  --  2.5,  h  -  2,  fl  -  0.5,  W'e  --  1,  zero  shear  rate  viscosity  - 10. 

Fig. 9 

Graph  of  g(z/h)  versus  z/h,  s  =  q  -  0. 

A  -  0.5.  C  -  2.5,  h  =  2,  fl  =  0.5,  We  -  1,  zero  shear  rate  viscosity^  10. 
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We  discuss  the  flow  of  BKZ  fluids  in  an  orthogonal  rheometer .  Some 
analytical  results  are  proved,  and  numerical  solutions  are  obtained  for  the 
Currie  model.  These  solutions  show  a  boundary  layer  behavior  at  high  Reynolds 
numbers  and  the  possibility  of  discontinuous  solutions  or  nonexistence  at 
high  Weissenberg  numbers. 
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